Bayesian spatial analysis of incomplete vaccination among children aged 12–23 months in Nigeria

High childhood disease prevalence and under-five mortality rates have been consistently reported in Nigeria. Vaccination is a cost-effective preventive strategy against childhood diseases. Therefore, this study aimed to identify the determinants of Incomplete Vaccination (IV) among children aged 12–23 months in Nigeria. This cross-sectional design study utilized the 2018 Nigeria Demographic and Health Survey (NDHS) dataset. A two-stage cluster sampling technique was used to select women of reproductive age who have children (n = 5475) aged 12–23 months. The outcome variable was IV of children against childhood diseases. Data were analyzed using Integrated Nested Laplace Approximation and Bayesian binary regression models (α0.05). Visualization of incomplete vaccination was produced using the ArcGIS software. Children’s mean age was 15.1 ± 3.2 months and the median number of vaccines received was four. Northern regions contributed largely to the IV. The likelihood of IV was lower among women aged 25–34 years (aOR = 0.67, 95% CI = 0.54–0.82, p < 0.05) and 35–49 years (aOR = 0.59, 95%CI = 0.46–0.77, p < 0.05) compared to younger women in the age group 15–24 years. An increasing level of education reduces the risk of odds of IV. Other predictors of IV were delivery at the health facility (aOR = 0.64, 95% CI = 053–0.76, p < 0.05), and media exposure (aOR = 0.63, 95%CI = 0.54–0.79, p < 0.05). Mothers’ characteristics explained most of the variability in the IV, relatively to smaller overall contributions from the community and state-level factors (p < 0.05). The level of IV against childhood diseases was high in Nigeria. However, disparities exist across the regions and other socioeconomic segments of the population. More efforts are required to improve vaccination sensitization programs and campaigns in Nigeria.

www.nature.com/scientificreports/(EPI), to cover preventable childhood diseases such as tuberculosis, poliomyelitis, diphtheria, pertussis, tetanus, hepatitis B, yellow fever, and measles.The EPI emphasizes that a child receives 1 dose of Bacille Calmette Guerin (BCG), 4 doses of oral polio-containing vaccine (PCV), 3 doses of diphtheria, pertussis, and tetanus (DPT) vaccine, and 1 dose of measles-containing vaccine (MCV).Additionally, the Polio Eradication Programme, yellow fever, Measles, and Meningitis Campaign programs were held recently in some states in Nigeria and this has increased the immunization levels in recent times.Although, Nigeria has made some progress in terms of reduction in under-five mortality which could be attributed to an increase in childhood vaccine uptake, about 4.3 million children still miss out on vaccinations yearly 3 .About 17.1 million children did not receive an initial dose of the DTP vaccine and an additional 5.6 million are partially vaccinated, and more than 60% of these children live in 10 countries including Nigeria 2 .The percentage of children aged 12-23 months who received all basic vaccinations increased from 25% in 2013 to 31% in 2018 5 .According to 2018 NDHS report, the prevalence of complete, incomplete and non-immunization of childhood vaccination in 2013 are respectively 25%, 54% and 21% while those reported in 2018 are respectively 31%, 49.9% and 19% 5 .These trends though indicate relative progress over the previous years.but are still far behind SDG3, for which the target is achieving more than 90% coverage of all basic vaccinations among children aged 12-23 months.
Previous studies [6][7][8] on immunization issues are mostly based on logistic regression without accounting for geographical clustering.Based on the efforts national and international agencies have put into ensuring national coverage, it is pertinent to assess the incompleteness of childhood vaccines uptake, identify geographical locations lagging, and understand the effect of factors on vaccination patterns across the geographical location, to enhance economic productivity.Therefore, this study took into consideration, mothers who reside in states where exposure to unmeasured similar social environment and health services access occur, so it is possible that there might exists spatial correlations across different communities, states and geopolitical regions.
The objectives of this study are to: examine the spatial patterns of vaccination coverage among children 12-23 months of age and identify the individual, community, and state-level factors associated with incomplete vaccination.The diverse cultural groups in Nigeria have implication for access to health care including childhood vaccination.The 1996 trial in Kano where five of the children given trovafloxacin and six who were given ceftriaxone died remains a barrier to childhood immunization uptake in the Northern part of Nigeria.Understanding the spatial variation in vaccination coverage in the context of states or administrative units covered is important for evaluating the performance of our vaccination programs.Thus, Bayesian regression models were employed to quantify the spatial risk of childhood vaccination as well as associations between incomplete childhood vaccination and a range of factors, and other fixed effects on childhood vaccination.
The conceptual framework for this study was based on Andersen's healthcare utilization and the Ecological Model as depicted in Fig. 1.Andersen Model was created to address disparities in healthcare access from the perspective of individual decisions made as a result of socioeconomic disparities 9 .The model posited that the usage of health services is determined by three dynamics: predisposing factors, enabling factors, and need.Predisposing factors can be characteristics such as race, age, and health beliefs.The enabling factors could be family support, access to health insurance, one's community, and so on.Need represents both perceived and actual need for health care services.This model has been used in the past for predicting healthcare utilization among population subgroups, and it takes into consideration of the minority groups, the rural dwellers, and the population group who may have challenges accessing healthcare.Similar to this model is the Social-Ecological Model which addresses health behaviors from the perspective that the observed health outcomes can be attributed to factors that may be analyzed at multiple levels which may go beyond individual characteristics.It considers the complex interplay between individual, relationship, community, and societal factors.This model suggests that interventions aimed at different levels may have a greater influence on health outcomes than interventions instituted at a single level 10 .
Figure 1.Conceptual framework adopted from Andersen's health care behavioral model.

Sampling technique
A two-stage cluster sampling method was used for data collection.In the first stage, 1400 clusters were selected with probability proportional to cluster size.The survey was successfully carried out in 1389 clusters across the country, a household listing procedure was carried out in all the selected clusters, and the resulting lists of households were used as a sampling frame where households were randomly selected for the interview.In the second stage, 30 households were selected in every cluster using systematic sampling, resulting in a total sample size of 5475.Detailed information about the sampling and data collection methods used during the survey is available at https:// dhspr ogram.com/ publi catio ns/ publi cation-fr359-dhs-final-repor ts.cfm.

Data source and extraction
The 2018 NDHS dataset, which is publicly available to all registered users, was requested from the Measure DHS (Demographic Health Survey) program website (https:// www.dhspr ogram.com).The available information about basic childhood vaccinations and their potential predictor variables were extracted.The geo-package of Nigeria's administrative areas at all levels of sub-division was used to identify hot spots of incomplete vaccination.

Outcome variable
The outcome variable was immunization status among children aged 12-23 months.According to the vaccination schedule in Nigeria, a child is considered fully immunized if he or she has received Bacille Calmette-Guerin (BCG) vaccination against tuberculosis and one dose of oral polio vaccine which are given at birth; three doses of DPT vaccine to prevent diphtheria, pertussis, and tetanus; another three doses of oral polio vaccine (OPV); and one dose of measles vaccine.These vaccines should be received during the first year of life 11 .Although, it is expected that within the first 12-23 months of life, all basic vaccinations should have been received by a child and this child is said to be completely vaccinated if he had received all of the immunizations recommended in the official immunization schedule 2 .To determine the immunization status, a new variable was generated for each of the vaccines and was regrouped as 1 for the uptake of a particular vaccine and 0 if a such vaccine was not received.Another variable was created to depict the total vaccine uptake for each child as follows: Lastly, any child that had received all the 9 basic childhood vaccines was classified as "1", incompletely immunized as "2" and not immunized as "3".

Explanatory variables
The basic explanatory variables are the individual and community-level characteristics.Individual-level variables included socio-demographic and economic characteristics such as the mother's age in years, mother educational level, mother's religious affiliation, mother's occupational status, place of delivery, the total number of children ever born, wanted the last child, exposure to media, birth order, wealth quintile, sex of the child, number of antenatal and postnatal care attendance within 2 months.Exposure to media, in this study, was created from the combination of questions on the frequency of mothers' access to information through newspapers/magazines, radio, or television.A set of variables including household assets, amenities and housing characteristics such as owner of television, refrigerator, type of flooring, source of drinking water, toilet facility, electricity, cooking gas and others were identified.The wealth index was constructed using principal component analysis (PCA).Each household in the dataset was assigned a wealth index score.These scores were then ranked and divided into quintiles and re-categorized into three wealth status; poor, average and rich.Community-level variables included the common characteristics of study subjects in an enumeration area such as place of residence,

Data analysis procedure
The data was weighted due to the complex sampling strategy used for data collection and to compensate for the unequal probability of selection between clusters.Descriptive statistics were used to describe the data.Highresolution maps of coverage of each vaccine and dropout rates between doses (for vaccines that exist in doses) were used to assess geographical variation in access to and acceptance of routine immunization.These maps were also integrated with high-resolution population maps to produce estimates of the numbers of undervaccinated children.SaTScan ArcGIS PRO 2.8 was used for spatial analysis and Stata 16.0 and R for non-spatial statistical analysis.SaTScan uses a circular window that moves systematically throughout the study area to identify statistically significant SaTScan clustering of areas with the same childhood basic vaccination coverage.Chi-square was used to examine factors that are associated with vaccination among children.At the level of multivariate analysis, binary logistic regression was used to identify factors influencing incomplete vaccination among children aged 12-23 months.In this study, the categories were dichotomized as 1 (if the child had completed the required vaccination) and 0 (if otherwise).A check of multicollinearity of the explanatory variables was carried out, and a condition number of 2.48 was obtained.This clearly indicated that the predicting factors have a low-level multicollinearity and are linearly independent.Thus, multilevel logistic regression models were generated and fitted to identify community and individual-level factors associated with childhood basic vaccination.The regression equation governing the use of multilevel modeling is given by Eq. (1): where; β 0j = γ 00 + γ 01 Z j + u 0j and β ij = γ 10 + γ 11 Z j + u 1j .
γ 00 and γ 01 are the intercept and slope to predict β 0j from Z j , u 0j is the residual error term in the equation for β 0j , γ 10 and γ 11 are the intercept and slope to predict β ij from Z j and u ij is the residual error term in the equation for β ij .Three models were fitted: the first model with individual-level variables (model I), the second model with individual and community-level variables (Model II) and the third model (Model III) with all the individual, community, and state-level variables.These models are given in Eqs.(2-4) Model comparison was carried out using Bayesian Information Criteria (BIC) and Akaike's Information Criterion (AIC).Finally, the model with the smallest value of the information criterion was selected as the best-fit model.Moreover, for measures of association (fixed effect), an adjusted odds ratio with 95% confidence intervals was used to determine statistical significance.For measures of variation (random effects), intra-class correlation coefficient (ICC), median odds ratio (MOR), and proportional change in variance (PCV) statistics were computed.ICC is a measure of within-cluster variation, the variation between individuals within the same cluster, and it was calculated using the formula: , where V A is the estimated variance in each model.The total variation attributed to individual or/and community level factors at each model was measured by the proportional change in variance (PCV) and was calculated as V A , where V A = variance of the initial model, and V B = variance of the model with more terms 12 .

The Bayesian Binary Regression model
The binary regression model 13 is used to explain the probability of a binary response variable as function of some covariates.The model is specified by: where y i = 1 if the response of interest is observed for the ith individual and zero otherwise, π i is the probability that the ith individual presents the response under investigation, β is the K vector of unknown parameters, X t i = x i1, x i2, . . ., x iK, the K vector of known covariates associated to the ith individual and F any transformation assuming values in (0, 1).The function F, for example, could be any arbitrary cumulative distribution function.The logistic, normal, and extreme values are the most useful F specifications, which are: (1) The link function defines the linear predictor as: For each of the three models stated above the linear predictor is respectively: . The model choice depends on the relation between the response variable and the covariates.
The likelihood function for data y = y 1 , y 2 , . . ., y n t is It is important to establish a joint prior distribution over the parameter space in order to proceed with the Bayesian analysis.Because the relationship between the data and the parameters is so intricate, this is extremely difficult to accomplish.The easiest way to circumvent this difficulty is to propose an informative prior, but with small precision, avoiding any complaint about the specification of subjective beliefs 14 .In this study, informative independent normal priors, with extremely small precisions, was set to the parameters.Therefore, Clearly, 9 is a complex function of the parameters and numerical methods are needed in order to obtain the marginal posterior distribution for each of the model parameters.Simulation based methods have proliferated in the last ten years or so yielding two popular approaches known as importance sampling 15 and Gibbs sampling 16,17 .

Ethical consideration
The survey was conducted under relevant guidelines and regulations.Permission to use the data was sought and granted by the data originator.However, ethical approval to conduct the study was obtained from the National Ethical Review Committee (NREC), and informed consent was also obtained from the respondents before the conduct of the interview.The respondents were assured of the anonymity of the information they provided.The possible identifier that could be used to track each respondent to the information they provided was removed from the data before use.

Results
The descriptive analysis of respondents by the vaccination status of children showed that a total of 5475 children were included in the study.1074 (19.6%) children did not receive any of the recommended vaccines, 3053 (55.8%), representing more than half of the children were partially vaccinated and 1348 (24.6%) of the children were only fully vaccinated.Figure 2 presents the percentage distribution of children who took each of the vaccinations under consideration.About three-quarters (72.9%) of the children received the second dose of the Polio vaccine and more than half (52.7%) of the children did not receive Polio 3-the fourth dose of the Polio vaccine.Also, 50.2% of the children did not receive the third dose of Diphtheria Tetanus Pertussis (DPT3).
Except for the sex of the child, Table 1 revealed that the vaccination status of the children was significantly associated with all of the explanatory variables.

Spatial prevalence of childhood vaccines uptake in Nigeria
Figure 3 disaggregated the spatial prevalence of childhood vaccination.No state was found to have a 0% prevalence but Kebbi, Sokoto, and Zamfara states in the Northern region had the least (˂5%) prevalence of complete vaccination.The highest prevalence above 50 was observed in Lagos and Anambra states.

Optimized hotspot analysis
The prevalence of incomplete vaccination among children aged 12-23 months in Nigeria is revealed in Fig. 4. Identifying clusters of various locations with a significantly high prevalence of incomplete vaccination is provided (7)

Random-intercept model of factors associated with childhood incomplete vaccination in Nigeria
The Multilevel Logistic model result is presented in this section.Model I (in Table 2) contains only the individual level variables, Model II contains an individual level and community level variables and Model III contains all levels' variables.The AIC and BIC of the three models show that the difference between them is quite negligible as they are very close.Even though Model I is the lowest in terms of AIC and BIC but Model III will be used for the interpretation since it captures more information.Mother's age and education, religion, place of delivery, media exposure, household wealth status, frequency of antenatal visits, postnatal care attendance within 2 months of delivery, community poverty level, state health facility, and population were all found to be significantly associated with childhood incomplete vaccination.Children whose mothers' ages fall between 35 and 39 years were found to have higher odds (OR = 1.66, 95% CI: 1.48, 1.92) of complete vaccination compared to children whose mothers fall in the age group 15-24 years.Women with primary and secondary/tertiary education were found to have higher odds [(OR = 1.15, 95% CI: 1.03, 1.50) and (OR = 1.51, 95% CI: 1.15, 2.40) respectively] of complete vaccination for their children compared to women who had no formal education.Children whose mothers practice none of Christianity or Islam religion were found to have lower odds of complete vaccination compared to those who practice Christianity (OR = 0.06, 95% CI: 0.01, 0.55).
Likewise, children who were delivered in the hospital/health center were found to have higher odds of (OR = 2.48, 95% CI: 1.80, 3.43) of complete vaccination relative to children who were delivered at home or elsewhere.Women who were exposed to media were also found to have higher odds of 40% (OR = 1.40, 95% CI: 1.27, www.nature.com/scientificreports/10.58) of complete vaccination for their children compared to women who were not exposed to media.Antenatal and Postnatal care attendance was also found to be significant.Children whose mothers went for 1-3, 4-7, and above 8 antenatal visits were found to have higher odds [(OR = 1.47, 95% CI: 1.12, 1.95), (OR = 1.47, 95% CI: 1.22, 2.47) and (OR = 2.85, 95% CI: 1.95, 3.81) respectively] of complete vaccination compared to children whose mothers did not go for any antenatal visit.Those who went for postnatal care attendance within two months of

Integrated nested Laplace approximation diagnostics on factors associated with vaccination status among children aged 12-23 months in Nigeria
Integrated nested Laplace approximation (INLA) offers an alternative to the computationally intensive Markov chain Monte Carlo (MCMC) methods 18,19 .INLA focuses on obtaining an approximation of the posterior marginal distributions of the model's parameters.To check whether the quality of the sample generated with the INLA algorithm was sufficient to provide an accurate approximation of the target distribution, the integrated nested Laplace approximation (INLA) diagnostics were employed.www.nature.com/scientificreports/Two diagnostics tools were used; the dispersion check and the distribution check.The classic dispersion measure is to sum the squared Pearson residuals and divide them by the number of observations minus the degrees of freedom.The dispersion check (D = 0.87) showed that there were no overdispersion and underdispersion, and no apparent anomalies in the chain, and the chain seemed to extensively explore the sample space many times.The distribution check takes the marginal distribution of the fitted values into account and returns the data frame with the empirical cumulative distribution of the raw data (ecdf) and the envelope of empirical cumulative distribution from data simulated from the model.The ratio of the ecdf of the raw data to the ecdf of the simulated data is near 1 (0.847).The model makes sense since the observed ecdf is within the envelope of the simulated ecdf.
Table 3 below presents the adjusted odds of factors associated with incomplete vaccination among children aged 12-23 months in Nigeria.Children of mothers aged 25-34 and 35-49 had lower odds of 33% and 41% (aOR = 0.67, 95% CI: 0.54, 0.82 and aOR = 0.59, 95% CI: 0.46, 0.77) to be incompletely vaccinated compared to those of mothers in the age group 15-24.There is also a lower odds of 28% (aOR = 0.72, 95%CI: 0.54, 0.95) of being incompletely vaccinated by children whose mothers had primary education compared to those children whose mothers had no formal education.
Similarly, the odds of childhood incomplete vaccination is 36% lower (aOR = 0.64, 95%CI: 0.53, 0.76) among children who were delivered in a hospital or other health center compared with children who were delivered at home or elsewhere.Children whose mothers were exposed to media (television, radio, and newspaper) had a lower odds of 37% (aOR = 0.63, 95%CI: 0.54, 0.79) to be incompletely vaccinated compared to the children whose mothers were not exposed to media.Household wealth status where a child belongs is a significant factor associated with incomplete vaccination.Children who belong to a household whose wealth status category is rich had lower odds of 42% (aOR = 0.58, 95%CI: 0.38, 0.89) of being incompletely vaccinated compared to children with poor household wealth status.
Additionally, there are lower odds of 63% (aOR = 0.37, 95%CI: 0.27, 0.52), 73% (aOR = 0.27, 95%CI: 0.20, 0.36), and 68% (aOR = 0.32, 95%CI: 0.23, 0.44) among children whose mothers went for 1-3, 4-7 and more than 8 antenatal visits respectively for incomplete vaccine uptake compared to children whose mothers did not go for any antenatal visit.Similarly, children whose mothers attended postnatal care within 2 months of delivery had lower odds of 35% (aOR = 0.65, 95%CI: 0.54, 0.78) to be incompletely vaccinated compared to those children whose mothers did not attend any postnatal care within 2 months of delivery.Likewise, children whose mothers reside where the state population is greater than 6 million had lower odds of 64% (aOR = 0.36, 95%CI: 0.15, 0.90) of incomplete vaccination compared to children whose mothers reside where the state population is less than 3 million.

Posterior prediction of vaccination status among children aged 12-23 months in Nigeria
The graph below (Fig. 5) is the posterior predictive check, which helps in inspecting the postulated or predicted values of vaccine uptake among children aged 12-23 months in Nigeria, against the observed values from the data.After using the observed data and predictors, we can update our belief about the unknown parameter (keeping all variables constant).The result predicts a continuous high number of children with no vaccination (over one-fourth of the population), then a sharp drop in the number of people that received one vaccination.Hereafter, a slight increase in the number of children taking 2 or 3 vaccinations (about 10%), accompanied by a steep decline in later vaccinations.This is contrary to the increase observed in the data after the 5th dose of vaccination.

Discussion
The prevalence of incomplete vaccination in Nigeria was high but strikingly higher in the Northern regions than in the regions in the South.This is due to relatively high level of non-formal education and low level of exposure to media of mothers in the north as compared to those in the southern parts of the country.This reveals that the mothers' educational level as a predisposing factor according to the Andersen's healthcare utilization model, influences children vaccines uptake behaviour.The findings of Adedokun et al. 20 and Afolabi et al. 11 corroborate with this assertion.The study revealed that one in every four children aged 12-23 months was not fully vaccinated in Nigeria.This finding is in agreement with the outcomes of the studies conducted by Sato and Eze et al., where similar level was observed 6,7,21 .The high level of incomplete vaccination observed in Nigeria could be possibly due to the inaccessibility of vaccination officers to some areas and localities where security issues were of serious challenge.With Andersen's healthcare utilization model, an enabling factor like healthcare providers is also found to impact childhood immunization rate in Nigeria.The perception of people about vaccination can be another possible reason for low vaccination coverage in Nigeria.When compared to children delivered at home or elsewhere, children delivered at health facilities receive all immunizations at birth; this is also reflected in an increased incidence of vaccination among women who attended antenatal care, and most importantly, postnatal care attendance within two months.Attendance in antenatal and postnatal clinics may be connected to a higher level of trust and confidence in the healthcare system.This finding is also in agreement with the study by Asuman et al. 22 who found that religion and household wealth are significant drivers and predictors of immunization and other health-seeking behavior.
The Bayesian output revealed that the states in the Northern part of Nigeria are mainly hot spots for incomplete vaccination.Although the distribution of incomplete vaccination spread across the country, it is worthwhile to note that childhood vaccination coverage is fair in the South-west.This study found that individual characteristics, community factors, and state factors all had a significant role in determining the variation in children's incomplete immunization status in Nigeria.Children of young mothers (15-24 years) are more likely to be under-immunized as compared to children of older women.This is in line with earlier studies 6,21,[23][24][25]  The study establishes the significance of education in predicting childhood immunization.The likelihood of having a child who is not completely vaccinated falls for a mother who just has elementary education or none at all compared to a child whose mother has secondary or higher education.It has been demonstrated that education significantly affects a mother's health-seeking behavior, including child vaccination 1,26,27 .According to comparable research, household poverty has an impact on children's vaccination rates 1,6,21 .A family's children are more likely to be under-immunized and unprotected against illnesses that can be prevented by vaccination the poorer the household grows.
The results of this study also demonstrate that pregnant women who use healthcare facilities help to raise the vaccination rates of their children.Children whose mothers received antenatal care had a decreased chance of not receiving the recommended vaccinations than children whose mothers did not get any prenatal care.This is also the result of earlier studies 27,28 .Attending prenatal care puts women in a better position to learn enough about routine childhood vaccination.
A mother also goes through processes during antenatal care that educates her to have a positive attitude toward using health care, not just for herself but also for her children.According to earlier studies, health facility delivery has a significant impact on baby immunization 1,6,21,29 .Incomplete vaccination is less common in children whose mothers gave birth in a hospital.Women who give birth in a health facility can learn about their child's immunization regimen and get their infant immunized at birth.
The ease of access mothers had to medical facilities was associated with childhood immunization.As previously reported, children of mothers who had trouble attending medical facilities are more likely to receive only partial vaccinations 8,21,23,30 .Children whose mothers live in socioeconomically poor regions and states are www.nature.com/scientificreports/more likely to be partially immunized.The findings of earlier studies 6,21,23 indicating the socioeconomic makeup of communities affected people's health-seeking behavior support this.

Limitations of the study
The study relied on secondary data.It is vulnerable to bottlenecks, omission mistakes, and manipulations because it is a secondary and retrospective cross-sectional survey (recall bias and no means of verification of all information supplied).Mothers' recall bias, which is often induced by the lack of vaccination cards, could lead to underestimating the study, albeit the NDHS methodology would mitigate the implications of these estimations.Vaccination coverage rates, on the other hand, are usually a measure of completing the specified vaccination schedule rather than of protection (the indicator does not reflect if the vaccines were given at the recommended ages and recommended intervals).Thus, the temporal distribution as well as the hierarchical nature of the number of vaccines taken and associated factors were not accounted for in the study.

Conclusions
This study considered the INLA random effects model to examine the factors associated with incomplete vaccination among children aged 12-23 months in Nigeria using the 2018 NDHS dataset.The performance of the model depended on the number of recommended vaccine uptake received by the children measured at different levels.A considerable number of children are not fully vaccinated against childhood infections in Nigeria.Individual and contextual factors significantly contributed to incomplete vaccination among children in the country.Therefore, the government and other stakeholders should pay particular attention to antenatal and postnatal attendance counseling programs for child-bearing mothers in all urban and rural parts of the country to overcome the problem of children's incomplete vaccination and/or vaccination dropout.Vaccines are a successful method for reducing child morbidity and mortality in nations all over the world, yet vaccination rates in Nigeria are extremely low.To boost the odds of uptake after birth, programs focused on encouraging and educating pregnant women, as well as young women, should be introduced.To meet the aforementioned recommendations, the government should enlist the help of non-governmental organizations and other relevant groups in promoting and maintaining vaccination uptake among children, particularly in low-uptake areas.

Figure 5 .
Figure 5. Plot of posterior predicts of the number of vaccinations among children aged 12-23 months in Nigeria (2018).

Table 1 .
Prevalence of vaccination status by socio-demographic factors in Nigeria.

Table 2 .
Factors associated with childhood incomplete vaccination in Nigeria.Significant values are in italics.
Vol.:(0123456789) Scientific Reports | (2024) 14:18297 | https://doi.org/10.1038/s41598-024-57345-y . This indicates that when a mother's age increases, the chance of her child receiving immunizations increases.The young mother's inexperience with child care may be to blame for this issue.Older women have faced the strains of caring for ailing children, with the ensuing effects on time and household finances.These mothers would support any effort to stop childhood illnesses.

Table 3 .
Bayesian model of factors associated with childhood incomplete vaccination in Nigeria.a Adjusted odds ratio significant at 95% confidence interval.